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Abstract 

Around 1950, Wigner introduced the idea of modelling physical re- 
ality with an ensemble of random matrices while studying the energy 
levels of heavy atomic nuclei. Since then, the field of random-matrix 
theory has grown tremendously, with applications ranging from fluc- 
tuations on the economic markets to complex atomic spectra. The 
purpose of this short article is to review several attempts to apply the 
basic concepts of random-matrix theory to the structure and radiative 
transitions of atoms and ions, using the random matrices originally 
introduced by Wigner in the framework of the gaussian orthogonal 
ensemble. Some intrinsic properties of complex-atom physics, which 
could be enlightened by random-matrix theory, are presented. 
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1 Introduction 

The structural and radiative properties of multi-charged ions are crucial for 
the studies of laboratory and astrophysical hot plasmas. Detailed numeri- 
cal calculations, however, are complicated and even sometimes impossible. 
Indeed, when the thermodynamical conditions are such that several shells 
are open, the total number of electronic configurations and electric-dipole 
lines can become huge. In that case the computation becomes numerically 
untractable. Complex atomic spectra pQ may be characterized by statistical 
(global) methods, particularly if one is interested only in general properties 
and regularities. The atomic spectrum as a distribution of energy levels 



or spectral lines can for instance be described by its statistical moments, 
which define its center of gravity, width, asymmetry (skewness), sharpness 
(kurtosis) and other mean characteristics [2]. A limited number of moments 
enables one to approximately describe the complex spectra with unresolved 
transition arrays [3], to reveal the statistical properties of spectra, to analyze 
the relative influence of various interactions, to study changes in atomic or 
isoelectronic sequences, to investigate configuration mixing, etc. The first 
problem is to find explicit expressions of such average characteristics [H [5] , 
without any detailed calculation of the spectrum. The second problem is the 
choice of the modelling function [6j [7] . Having explicit formulas for a number 
of first moments enables one to approximately restore the envelope of the 
radiation spectrum. Studies of variation of these statistical characteristics 
along isoelectronic sequences give a wealth of information on intra-atomic 
interactions [8]. 

However, such a statistical approach does not provide any information 
about the internal correlation laws of a transition array. It gives no indica- 
tion about the possible occurence of chaotic phenomena. 

The random- matrix theory (RMT) [9l fl0l[TT| ll2j is a stochastic approach 
which relies on the assumption that the Hamiltonian matrix represents an 
ensemble of real finite matrices characterized by the values of their elements, 
which are defined by their differential probability to occur. The matrix 
probability distribution is assumed to consist of a product of distributions 
for the individual matrix elements (i.e. the matrix elements are statistically 
independent). Wigner's gaussian orthogonal ensemble (GOE) consists of 
real symmetric matrices where the distribution of the diagonal elements is 
given by a Gaussian with a zero mean and a variance twice that of the off- 
diagonal elements. Obviously, the matrix elements of the real system are 
not distributed independently. However, the validity of the approximation 
can be tested on the predictions made for the eigenvalue spectrum. RMT 
predicts the way adjacent levels are correlated (statistical "repulsion" of 
levels) . 

The successful application of RMT to the study of the distribution of 
energy levels of atoms and ions is recalled in section O and its extension to 
the calculation of electric-dipole (El) lines in section [3l An interesting al- 
ternative, the Two-Body Random Matrix Model is also explained in section 
El In section HI several unexplained properties of complex atomic spectra 
are presented, within an attempt to relate them to RMT. Section [5] is the 
conclusion. 

2 Energy levels 

In the atomic central-field model, the Hamiltonian is diagonalized in elec- 
tronic configurations. The distribution of the matrix elements is so complex 
that it can be considered as random. In the RMT, all the elements of a 
real symmetric matrix in Wigner's GOE are independent random variables 
obeying a Gaussian distribution. The eigenvalues of a random matrix are 
highly correlated; for instance, two adjacent eigenvalues "repel" each other. 
The distribution of the eigenvalues themselves tends to become Gaussian 
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when the number of interacting particles increases. In contrast, numbers 
chosen at random are not correlated. More precisely, if the level spacings 
are considered as a set of numbers Xj+i = Xi + s with x% = such that s is 
a random variable [13], their distribution obeys Poisson's law: 



With random correlations of the matrix elements s of H , the level spac- 
ings follow the Wigner surmise 



which holds in the GOE. Actually, the GOE represents a special case of 
a large family of ensembles which all have the same spectral behaviour. 

Rosenzweig and Porter [14] plotted the empirical distribution of nearest- 
neighbor spacings for the odd-parity levels of neutral Hf (Z=72). The 
Wigner distribution (see equation ([2])) was found to give an excellent qual- 
itative agreement with the empirical distribution. Later, Flambaum et al. 
[15] studied the spacing distributions of the J 71 = 4 + levels of Ce (Z=58). 
They found that the experimental level spacings were in good agreement 
with the Wigner distribution. O'Sullivan et al. |16j pointed out that taking 
configuration interaction into account is crucial in order to obtain a Wigner 
distribution for the level spacings. The single-electron states mix together 
and produce a spectrum of levels analogous to the one in compound nuclei. 
In fact, if the ratio of the mean off-diagonal Coulomb interaction parameter 
by the mean level spacing is large, the single-particle model no longer holds 
and the states highly mix so that only total angular momentum J and parity 
7r are good quantum numbers. 

Studying the distribution of the adjacent-level spacings enables one to 
check how far the level density lies from the Poisson and Wigner (derivative 
of a Gaussian) limits [T71 [TS] . In the GOE, a \ 2 distribution can be derived 
from the Porter-Thomas law [191 120] for the strength of the lines (see section 



Bohigas, Giannoni and Schmidt [21] conjectured that the spectrum of 
a chaotic system is strongly related to the spectrum of a random matrix. 
More precisely, they found that the level fluctuations of the quantum Sinai's 
billiard are consistent with the predictions of the GOE, which was a confir- 
mation that RMT is strongly linked to quantum chaos and that the fluctu- 
ation laws are universal. In fact, for classical systems, non-integrable pro- 
cesses follow Wigner distribution, while the properties of integrable systems 
obey Poissonian statistics. The extrapolation from classical chaotic systems 
to quantum ones was discussed for instance by Connerade et al. [22j [23] . 
Yukawa recovered the GOE predictions without introducing a priori ensem- 
bles of random matrices. By means of the usual statistical laws of classical 
many-body systems, he proposed a level distribution function which makes 
a transition from the Poisson type to the Gaussian type depending on the 
value of a single parameter characteristic of the Hamiltonian. 
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3 Electric-dipole lines 



While the position of the absorption lines depends on the eigenvalues, their 
strength is mostly governed by the eigenvectors. In an atomic-structure 
code, most of the time is consumed in diagonalizing the Hamiltonian matri- 
ces and constructing the dipole transition matrix. Most of the calculational 
effort resides in constructing the off-diagonal matrix elements of the blocked 
Hamiltonian matrices. This is easily understood as these matrix elements 
are formed not only from one-body operators (as the transition matrix ele- 
ments) but also from two-body Coulomb interaction operators. 

3.1 The RMT as a tool for simplifying the exact calculation 

In their work [25j, Wilson et al. calculated the diagonal terms of the Hamil- 
tonian matrix in a pure coupling using Cowan's atomic-structure code |26j 
and populated off-diagonal elements statistically beyond the GOE according 
to a bi-Gaussian distribution function. They observed a disproportionately 
large number of off-diagonal elements of small amplitude and noticed that 
the larger matrix elements connect basis states where the parent shell is also 
of common spectroscopic term. A detailed analysis showed that the larger 
elements are distributed according to the GOE predictions, while the smaller 
are distributed like a Gaussian with a much smaller width (bi-Gaussian sur- 



3.2 The two-body random interaction model 

The RMT does not describe many important properties of realistic many- 
body systems, its predictions being limited to level statistics, localization 
properties of the eigenstates, etc. The Two-Body Random Interaction Model 
(TBRIM) [271 EHl [291 SO] was introduced as an intermediate approach based 
on a simple mathematical model with random interactions which, however, 
takes into account the most important features of many-body systems: or- 
bitals, two-body interactions and the Pauli exclusion principle. It consists 
of a system of fermions (n fermions distributed in m orbitals with energies 
€i distributed between and 1, see Refs. [3TJ [32], [331 G2] ) ■ Then, all possible 

(k) 

states with n fermions are built, the state k being characterized by n\ = 
or 1, with i= 1, ■ • ■ , m. The Hamiltonian matrix is such that 



and Vki is distributed randomly in the interval [— Vo, Vq\ and is set equal 
to if k and / states differ by more than two orbitals. 5ki is Kronecker's 
symbol. The density of states is found to be a Gaussian: 
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centered in 
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1 N 

E c = ^Y. R kk (5) 
k=l 



and having a variance 

N -, N 



AT ^- — ' c /V 

fe=l fc=l 



where 



AT 

<A£ fc 2 > = £ H lv (7) 
The number iV is the size of the two-body random interaction matrix: 

V ra / n\(m — n)\ 
and the sparsity r of the matrix is given by [35j: 

r = — [1 + n(m — n) + n(n — l)(m — n)(m — n — l)/4]. (9) 

The diagonalization of the Hamiltonian matrix leads to the following 
expression of the wavefunctions: 

N 

IV> w ) = £cf \4>j), (io) 

where the coefficients C- , defined by 

are the components of the eigenstates in terms of the \4>j) basis. The 
Ej = Hjj expectation value represents the energy of the many-electron atom 
state \4>j) and 

N N 

7i = (A4)=£|cf| (Ej-EW) = £ H%, (12) 

i=l P=l,p¥=j 

where E^> is the energy of state \ifj^). The configuration interaction 
leads to the "straightening" of the spectrum. This effect is a direct result of 
the level repulsion. It strongly manifests in the parts of the spectra where 
the average level spacing is much smaller than the typical matrix element iTj, 
(i j). In complex atomic spectra, if a given interaction mixes efficiently 
the basis "natural" states, the states of the system are composite states with 
the fractional occupation numbers for the orbitals 

n s = (^\ata s \^) = J2\cl l) \ (faK a s \cf> k } , (13) 

k=l 
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where a + and a represent the second-quantization creation and annihila- 
tion operators respectively. Angom et al. [36J showed, through a statistical 
analysis of Sm I eigenvalues and eigenfunctions for J = + and J = 4 + , that 
atoms with a sufficiently large number of active electrons exhibit TBRIM 
characteristics. The density of states predicted by the TBRIM is Gaussian 
and a Gaussian-like density of states was also observed in the Sm IX spectra 
by O'Sullivan et al. [16J. The fact that some quantities can be obtained 
without diagonalization of the Hamiltonian matrix can be useful for the 
calculation of transition rates [Ml E7j . 

4 Intrinsic properties of spectra 

In this section we mention a few interesting properties of the transition 
arrays of El lines which are not fully understood yet and for which the 
RMT could bring some explanations. 

4.1 Scars of symmetries 

In transition arrays of complex atomic spectra, some well-defined subsets of 
lines bear the signatures, called scars, of their vanishing in one or several 
pure couplings. Symmetries are essential in the approach to n— particle 
systems. In atomic and molecular physics, some current studies deal with 
complex low-symmetry situations. For example, the high symmetry of the 
hydrogen atom, when broken by an external field, can lead to quantum 
chaos. However, it was shown by Delande and Gay [38J that, for that system, 
symmetry effects can also be detected, as scars in the joint distribution of 
line wave numbers and oscillator strengths. In their paper, Bauche and 
Bauche-Arnoult [39] presented some scars of symmetries in the distribution 
of line intensities in complex atomic spectra. They found that in a particular 
type of transition array, one can define subsets of lines whose lines are, on 
the average, definitely weaker or stronger than the others. The fact that 
they are weaker is trivial when the physical coupling is close to an extreme 
coupling where they vanish. This was shown numerically, by means of line- 
by-line calculations, in the cases of two particular transition arrays d 4 — > d 3 p 
and d 8 — > d 7 p. The authors pointed out that such a phenomenon can persist 
very far from that coupling, which was rather unexpected. Such scars are 
the traces of the mathematical structure of the intensity distribution, which 
remains observable in extremely complex situations. The corresponding 
scars are different from the ones exhibited by Delande and Gay [38], since 
they are not produced by an external field; they are observed in the intensity 
distribution without reference to the line wave numbers and they occur in 
many-electron atomic systems. 

4.2 Fractal nature of atomic spectra 

Learner [40j found what he called an "unexpected law" related to the number 
of weak lines in a spectrum. The logarithm log(A^) of the number of lines 
whose intensities lie between 2 k Io and 2 fc+1 /o is a decreasing linear function 
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of k. In the example of Learner (spectrum of Fe I), the value of Iq is chosen 
in such a way that this law holds for 1 < k < 9 (9 octaves) when about 1500 
lines within the wavelength range 290 nm < A < 550 nm are considered. If 
p is the slope, one has = Nq x 10~ fcp . The number of lines is divided 
by 10 p when the size of the interval is multiplied by two. A characteristic 
dimension d is therefore given by 

— = 2 d , (14) 

IQP ' v 1 

i.e. 

d = — - — — 7— r- (15) 

logio(2) 1 ' 

Learner noticed that p ~ log 10 (2)/2. In that case, d ~ —1/2 and it will 
be shown in a future work that the Hausdorff fractal dimension associated 
to Learners's rule is close to 1/3. The link between fractal behavior and the 
RMT is not obvious, but definitely exists [41]. An excess of weak lines is 
also observed in that case, as in the approach of Ref. [25]. Bauche et al. 
related this excess to the coupling effects [42] when the physical coupling is 
close to a pure coupling, and to scars of symmetries in the opposite case. 



4.3 The law of anomalous numbers 

The observation by Newcomb [43] that the first pages of logarithm books 
were more used than the last ones, led to the conjecture that the signifi- 
cant digits of many sets of naturally occuring data are not equi-probably 
distributed, but in a way that favors smaller significant digits. For instance, 
the first significant digit (i.e. the first digit which is non zero) will be 6 more 
frequently than 7 and the first three significant digits will be 439 more often 
than 462. Benford provided a probability distribution function for signifi- 
cant digits, which states that the probability that the first significant digit 
d\ is equal to k is given by [33] : 

P(d 1 =k) = log 10 (l + ^) • (16) 

It was found recently (see Ref. [H]) that the distribution of lines in a 
given transition array follows very well Benford's logarithmic law of signifi- 
cant digits. This indicates that the distribution of digits reflects the symme- 
try due to the selection rules. If transitions were governed by uncorrelated 
random processes, each digit would be equi-probable. 

Benford's law is still not fully understood mathematically. However, it 
applies if the system is governed by multiplicative processes [46]. Indeed, a 
random multiplicative process corresponds to an additive process in a loga- 
rithmic space. As we have seen before, in Wigner's RMT, the Hamiltonian is 
defined in the GOE by an ensemble of real symmetric matrices for which the 
probability distribution is a product of the distributions for the individual 
matrix elements H^i, considered as stochastic variables, and the variance 
of the distribution for the diagonal elements is twice the one for the off- 
diagonal elements. The line strength S is proportional to |(i|I)|j)| 2 , where 
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D is the dipole operator and \i) and \ j) are eigenvectors of the Hamiltonian. 
Therefore, the line strengths involve quantities such as products of H^i and 
one has 



log S = log C + log S', (17) 

where ( is a stochastic variable and S and S' two different strengths. 
The central-limit theorem states that the probability distribution that the 
value of the n th strength is S will be Gaussian with a variance oc n 1 / 2 . In 
the infinite limit, the distribution will approach a constant K, so that one 
has 

dS 



J P(logS)d(logS) = K J — . (18) 



The probability P that the first significant digit d\ of S is k in base 10 
is given by 

f k+1 dS I f l0 dS , / U 
m=t) = X -8 Ik -S =] ° g "'( 1 + kl (19) 

which is exactly equation (I16p . This argument relies on the idea that 
fluctuations are governed by multiplicative processes involving a stochastic 
variable. The matrix elements of the Hamiltonian are correlated stochastic 
variables and the product of such variables leads to Benford's logarithmic 
distribution of digits. 

Another argument can be invoked to understand the validity of Benford's 
law for the analysis of complex atomic spectra. As we have seen in section 
[2l Porter and Thomas have shown that the amplitudes of the lines between 
all the levels of two random matrices obey a Gaussian distribution [THJ [20], 
which is 



D{s) = vmm ew 



s 



2(5) 



(20) 



where L and (S) are respectively the number of lines and the average 
value of the line strength S (the strength is the square of the amplitude) in a 
(J, J') set. When numbers are taken from an exponential distribution, they 
automatically obey Benford's law. Therefore, when the exponential term 
dominates in D(S), which is often the case except close to the origin (i.e. 
for very weak lines), Benford's law applies. 

The RMT contains approximate symmetries, which are not sufficient 
to describe the vicinity of Russell- Saunders (LS) and jj couplings [39J 147] . 
However, we found that even when the RMT is expected to be inappropriate, 
i.e. close to a pure coupling, the line strengths still fit Benford's law. This is 
due to the fact that even if the number of weak emerging lines is important 
(when the term -7= dominates in the Porter-Thomas law (|2"0"])). the number 
of decades is sufficiently large so that the weight of those weak lines is 
negligible in the statistical occurence of digits. 

Since Benford's law can be explained in terms of a dynamics governed 
by multiplicative stochastic processes, the RMT is probably an interesting 
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pathway for the calculation of large atomic-dipole transition arrays and Ben- 
ford's law can help clarifying the existence of different classes of stochastic 
Gaussian variables. 

5 Conclusion 

Random-matrix theory offers a new kind of statistical description of disor- 
dered systems. It is concerned with the statistical properties of eigenvalue 
sequences. Since its inception in the field of nuclear physics, this descrip- 
tion was found to apply to the statistical description of the eigenvalues of a 
quantum chaotic system. Although the level sequences that arise from these 
different physical systems follow different overall distributions, the fluctua- 
tion properties that they exhibit seem to be universal. This universality is 
an important property of RMT but as with any statistical description, de- 
tailed information of a system is lost. The strength and relevance of RMT lie 
in the fact that the fluctuation properties of a system's eigenvalue sequence 
are independent of its detailed structure. Rather, the fluctuation properties 
of such a sequence depend solely on the system's overall symmetries (which 
are called "couplings" in atomic spectroscopy). RMT can therefore be used 
to get insight into the symmetry properties of a system by studying the sta- 
tistical properties of its corresponding eigenvalue sequence. Since the fifties, 
RMT has grown tremendously in scope, and has even become a field of 
study in its own right. With this growth, the mathematical tools developed 
allow ensembles of random matrices to model the structural and radiative 
properties of complex atoms encountered in multicharged-ion plasmas. 
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